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00 . To predict the X-ray observables associated with the diffuse baryons in clusters of 

^ \ galaxies, we develop a new physical model for such a hot intra-cluster plasma. Our 

• framework is provided by the hierarchical clustering cosmogony for the dark matter, 

and by the standard FRW or Lemaitre cosmologies constrained by cosmic ages. 
O I As to the plasma dynamics and thermodynamics we propose a semi-analytical ap- 

\ proach based on punctuated equilibria. This comprises the following blocks that we 

■ compute in detail: Monte Carlo "merging histories" to describe the dynamics of dark 

matter condensations on scales of order 1 — 10 Mpc, and the associated evolution of 
^ \ the gravitational potential wells; the central hydrostatic disposition for the ICP, reset 

^1 to a new equilibrium after each merging episode; conditions of shock, or of closely adi- 

abatic compression at the boundary with the external gas, preheated by stellar energy 
feedbacks. Shocks of substantial strength are shown to prevail at the outskirts of rich 
clusters in a universe with decelerated expansion. 

From our model we predict the L — T relation, consistent with the data as for shape 
and scatter. This we combine with the mass distribution provided by the canonical 
hierarchical clustering; the initial perturbation spectra are dominated by Cold Dark 
Matter but include enough baryons to account for the high abundance sampled by the 
X-ray clusters, and are COBE-normalized. Thus we predict the z-resolved luminosity 
functions, with the associated source counts and redshift distributions. We predict also 
the complementary contribution by the unresolved groups and clusters to the soft X-ray 
background. 

These results are compared with two recent surveys from ROSAT; one defines the 
local luminosity function over nearly three decades of L, and the other shows little or 
no evolution out to z ~ 0.8. 
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Our results confirm that the critical cosmology coupled with Standard CDM is ruled 
out by its overproduction of local clusters. On account of underproduction, instead, we 
rule out open cosmologies (the cheapest way to solve the baryonic crisis and to freeze 
evolution), except for a narrow range around flo — 0.5; even there, we find the con- 
sistency with the full data base to be hardly marginal. For the CDM cosmogony with 
flo = 0.3 but in fiat geometry, we obtain acceptable fits. For the tilted CDM pertur- 
bation spectrum with high baryonic content in the critical universe, we obtain marginal 
consistency. The cosmogonical/cosmological sectors of the cluster history are indepen- 
dently testable by means of a lower bound to the evolved temperature distribution, as 
can be measured with SAX and XMM out to moderate z. 

Finally, we discuss the effective limitations of X-ray clusters and groups as cosmo- 
logical signposts, and their brighter prospects toward the astrophysics of the ICP and 
the cosmogony of large, high-contrast structures. 

Subject headings: galaxies: clustering - galaxies: intergalactic medium - galaxies: X-rays 
- hydrodynamics 
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1. Introduction 

Groups and clusters of galaxies constitute cosmic 
structures siifficicntly close to eqiiilibrium and with 
sufficient density contrast (5 » 2 10^ inside the virial 
radius R) as to yield definite observables, and possibly 
to offer reliable signposts for cosmology. 

They stand out to substantial depths of space- 
time not only in the optical band, but even more 
in X-rays. This is because their gravitational po- 
tential wells, shaped by a dominant dark mass M, 
contain not only baryons condensed into stars but 
also a larger amount of diffuse baryons. The latter, 
with densities n ~ 10""^ cm~^ and virial temperatures 
kT 1^ GM niH/'i-OR ^ 5 keV in rich clusters, sat- 
isfy the plasma condition kT/e^ n}-!^ ^ 1 exceedingly 
well; in fact, they do by a factor 10^^, vastly larger 
than for the baryons inside the stars. Such a hot 
intra-cluster plasma (ICP) emits powerhil X-ray lu- 
minosities L (X -n? R\ T^/^ ~ 10'*^ erg/s by optically 
thin thermal bremsstrahlung from central regions of 
overall radius Rx ^ 1 ATpc. 

The temperature directly probes the height of the 
potential well, with the baryons in the role of mere 
tracers; on the other hand, the luminosity with its 
strong dependence on density reliably probes the 
baryonic content. Statistically, a definite L — T corre- 
lation is observed (albeit with considerable scatter), 
and this provides the crucial link to relate the X-ray 
luminosity functions with the statistics of the dark 
mass M or with that of the corresponding T. 

But groups and clusters are intrinsically complex 
systems. To begin with, their dynamical history is 
marked by extensive, repeated merging of clumps, 
both in the form of nearly isotropic accretion of small 
units, and in the form of a few large, anisotropic co- 
alescence events. This is shown in detail by all N- 
body simulations (see, e.g., Schindler & Mueller 1993; 
Tormen 1997; Roettiger, Stone & Mushotzky 1997), 
and is increasingly indicated by the data (see, e.g., 
Zabludoff & Zaritsky 1995; Henriksen & Markevitch 
1996; Jones et al. 1997). 

The timing of such a dynamical evolution (and 
specifically the present merging rate) is set by the 
cosmological framework. Adopting the homogeneous 
isotropic FRW or Lemaitre cosmologies, the expan- 
sion is parametrized by the Hubble constant h (in 
units of 100 km/s Mpc), and by the density parame- 
ter (to which the baryons contribute the fraction 
Ob); a possible additional contribution is given by Qx 



associated to the cosmological constant. 

Given the cosmological framework, the cosmogony 
(i.e., the process of structure formation) is treated 
in terms of the standard hierarchical clustering sce- 
nario (see Peebles 1993), where all structures form 
by gravitational instability of initial density pertur- 
bations S in the dark matter (hereafter DM). Impor- 
tant parameters are the shape of their power spectrum 
{\Sk\'^) oc fc"" T^(fc) at the recombination, and the nor- 
malization measured at large scales by COBE/DMR 
(Gorsky et al. 1996). The transmission function 
T{k) depends on the specific model assumed for the 
DM; given this, the amplitude as (i-c., the normal- 
ization extrapolated down to the relevant scale 8h~^ 
Mpc) and the effective spectral slope rig depend on 
rio and h. These parameters also affect the growth 
factor D{z) of the perturbations, which enters the 
actual predictions for the mass distributions of the 
condensed clusters or groups (see Appendix A). 

Further complexity is added by the physics of the 
IGP, and this constitutes our main aim here. To now, 
the ratio of the ICP to DM, and specifically its central 
density Uq and its effective radius arc poorly under- 
stood. But the former is especially important, as the 
factor rip enters the luminosities and amplifies the ob- 
served variance. 

Observational information on the ICP and on the 
underlying DM dynamics is provided by the local 
L — T correlation. At higher z, further statistical 
information is provided by the evolution of the X-ray 
luminosity function N{L,z), or by its integrals like 
the source counts or the z-distributions. The prob- 
lem is that the predictions of observables involve not 
only the ICP physics and the cosmogony (with their 
intrinsic variances) but also cosmology again (with 
its uncertainties), and so the various aspects are not 
easily disentangled. 

The sharpest result obtained so far rules out the 
attractively simple assumption by Kaiser (1986) that 
the ICP amount be proportional to the DM's from 
groups to clusters at all z and M. A large number 
of subsequent works (of which we cite here the recent 
Mathiesen & Evrard 1997, Kitayama & Suto 1997, 
and Borgani et al. 1998) dealt with the combinations 
of these three sectors, namely, cosmology, cosmogony 
and ICP physics. But most of these papers, following 
the start by Cavaliere & Colafrancesco (1988), ap- 
proached the problem by parameterizing the depen- 
dences of the ICP/DM ratio on M and z, e.g., in the 
form L oc (1 -|- zy. While the parameters p and s 
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arc constrained to some extent by the local L — T cor- 
relation, nevertheless there still remains a substantial 
degeneracy (Oukbir, Bartlett & Blanchard 1997) be- 
tween ICP physics and cosmology /cosmogony. That 
is to say, different combinations of ft, Qo, ^A, ^b, cts, 
Up, s, p provide close fits to the observables. Con- 
versely, when the ICP parameters are varied, differ- 
ent cosmologies appear to be preferred by the data; 
the trade off particularly concerns Qg and s, which 
directly govern the global and the ICP evolution, re- 
spectively. 

To go beyond such degeneracies a physical model is 
needed for the ICP, including the above complexities. 
Such a model must include: the histories of DM halos 
with their hierarchical merging events; the infall of 
the gas with the ensuing compression and shocks; the 
disposition of the ICP in the potential wells; its condi- 
tions at the boundary with the surrounding environ- 
ment, which is modulated by the large scale structures 
and by stellar preheating. A model accounting for all 
the above is missing so far. We stress that the simula- 
tions using advanced hydrodynamics coupled with N- 
body codes hardly reach at present enough dynamic 
range (as discussed by Bryan fc Norman 1997) to de- 
scribe DM and ICP over the full range from ~ 50 
Mpc associated with the large scale structures (which 
guide the ongoing mergers of DM halos) , to the inner 
50 kpc where the ICP yields a considerable contribu- 
tion to L. On the other hand, such high-resolution 
simulations do not include yet the stellar preheating 
with its crucial effects on the LT relation. 

This motivates us to develop here a semi-analytical 
model which includes, though in a simplified form, the 
features listed above. We describe the cluster evolu- 
tion as a sequence of punctuated equilibria (PE) ; that 
is to say, a sequence of hierarchical merging episodes 
of the DM halos (computed with Monte Carlo sim- 
ulations), associated in the ICP to shocks of various 
strengths (depending on the mass ratio of the merg- 
ing clumps), which provide the boundary conditions 
for the ICP to re-adjust to a new hydrostatic equilib- 
rium. We show that our PE model predicts density 
and temperature profiles and the L — T relation for 
clusters and groups consistent with the recent data. 

We then use our PE to predict the counts of 
resolved sources N{> F) for faint fluxes down to 
F > 2 10""'^^ erg/s cm^, now accurately measured by 
Rosati et al. (1997). We predict also the complemen- 
tary observable constituted by the contribution of the 
unresolved groups and clusters to the soft X-ray back- 



ground (XRB). We predict the z-distributions, and 
finally the full z-resolved luminosity functions to be 
compared with recent and forthcoming surveys. 

The paper is organized as follows. In §2 we present 
and discuss our approach to the ICP astrophysics. In 
§3 we give the X-ray observables in the form suited to 
our hierarchical clustering computations; details are 
supplied in the Appendixes A, B and C. In §4 we 
present the results from our approach. The final §5 is 
devoted to discussions and conclusions. 

2. The punctuated equilibria for the ICP 

The X-ray luminosity of a cluster with radial tem- 
perature profile T(r) and density profile n(r) is given 
by 

Loc J d'^rn^(r)T^'^{r) ; (1) 

the integration is over the emitting volume r^ < R^. 
Expressing R from T oc M/R oc pR^, eq. (1) recast 
into the form 

ioc [n(r)/ni]2T2pi/2 , (2) 

where p is the internal DM density, ni is the den- 
sity just exterior to the cluster boundary, and the 
bar indicates the integration over the cluster volume 
normalized to i?'^. Note that eq. (2) applies in the 
isothermal case; the corresponding expression for a 
polytropic ICP is given in Appendix B. 

The simplest approach to the ICP state is that 
adopted in the self-similar model (Kaiser 1986) where 
n (X p is assumed, independently of T; then L oc 
obtains from eq. (2). The result conflicts with the 
observed correlation for rich clusters which is close 
to L (X (Edge & Stewart 1991; Mushotzky 1994; 
Tsuru et al. 1996; Mushotzky & Scharf 1997). Also, 
when combined with the standard hierarchical cos- 
mogony, the assumption yields unacceptable fits to 
the local luminosity function (see, e.g., Kitayama & 
Suto 1997). Finally, it would predict the clusters at 
higher z to comprise not only denser DM, but also 
equally denser ICP; since bremsstrahlung depends on 
n^, this would imply a strong positive evolution of 
N(L,z), which is certainly not observed. Rather, 
the analysis of the ROSAT Brightest Cluster Sample 
(Ebeling et al. 1997) and of the ROSAT Deep Clus- 
ter Survey by Rosati et al. (1997), extending and 
strengthening the data by Collins et al. (1997) and 
by Nichol et al. (1997), indicate no evolution out to 
z w 0.8 for L < 3 10^^ erg/s; earUer surveys (EMSS, 
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Henry c't al. 1992) siiggc'stcd oven a (marginal) nega- 
tive evolution ol the bright clusters. 

So to derive the true scaling of L with T we need 
a closer analysis of the IGP disposition relative to the 
DM. We propose a new approach based upon two cor- 
nerstones: the profiles n(r) and T(r) are given by the 
hydrostatic equilibrium with the gradually changing 
gravitational potential; their normalizations, and so 
the central density, are set by the conditions at the 
boundary with the external medium. 

The equilibrium profile may be effectively repre- 
sented with a polytropic relation starting from the 
cluster boundary at r2, say: 



n{r)/ni = g{T) 



T{r) 



1/(7-1) 



(3) 



Here T2 is the temperature just interior to the bound- 
ary; we conveniently use g{T) = n2/"i for the ratio 
of the interior to the exterior density, to include any 
shock discontinuity at the boundary. The appropriate 
values for 7 will be discussed in §2.1 and §2.4. 

Actually, the ICP is reset to a new equilibrium after 
each episode of accretion or merging of further mass. 
In our PE approach, the history of such episodes is 
followed in the framework of the hierarchical cluster- 
ing by Monte Carlo simulations, as explained below. 

Two relevant limiting forms of eq. (3) are consti- 
tuted by by the "shock" model of CMT97, and by the 
"adiabatic" models of Kaiser (1991) and of Evrard 
& Henry (1991). In both the outer gas is expected 
to be preheated at Ti ^ 1 keV (Ciotti et al. 1991; 
David et al. 1993, 1995; Renzini 1997) hy z <^ 2, 
due to feedback energy inputs following star forma- 
tion and evolution. Preheating temperatures Ti ^0.1 
keV also would prevent the cooling catastrophe from 
occurring, see White & Rees (1978); Blanchard, Vails 
Gabaud & Mamon (1992). In point of fact, Henriksen 
& White (1996) find from X-rays evidence for diffuse 
cool gas at 0.5 — 1 kcV in the outer regions of a num- 
ber of clusters. In the present context, preheating 
will inhibit the attainment of the universal baryonic 
density in gravitational wells with virial temperatures 
comparable to Ti. 

These limiting models differ in their treatment of 
the boundary conditions and of the merging histories. 

2.1. Shocks and hydrostatic equilibrium 

The key boundary condition is provided by the dy- 
namic stress balance P2 = Pi+niH n\ vf, relating the 



exterior and interior pressures P2 and Pi to the inflow 
velocity vi driven by the gravitational potential at the 
boundary. We expect the inflowing gas to become su- 
personic in the vicinity of R, when mjj > 2kTi. 
In fact, many hydrodynamical simulations of loose 
gas accretion into a cluster (from Perrenod 1980 to 
Takizawa & Mineshige 1997) show shocks to form, 
to convert most of the bulk energy into thermal en- 
ergy, and to expand slowly remaining close to the 
virial radius for some dynamical times. So we take 
r2 « i? (which follows the structure growth, since 
R oc M^/^), and focus on nearly static conditions in- 
side, with u| << vf. 

The post-shock state is set by conservations across 
the shock not only of the stresses, but also of mass 
and energy, as described by the Rankine-Hugoniot 
conditions (see Appendix B). These provide at the 
boundary the temperature jump T2/T1, and the cor- 
responding density jump g which reads 



9{ 



Ti 
T2 



1/2 



(4) 



for a plasma with three degrees of freedom. Eq. 
(4) includes both weak and strong shocks. For weak 
shocks with T2 « Ti (appropriate for small groups 
accreting preheated gas, or for rich clusters accreting 
comparable clumps), this converges to the truly adi- 
abatic relationship n2/ni — {T2/Ti)^/'^ up to second 
order inclusive, see Landau & Lifshitz (1959). On the 
other hand, it is shown in Appendix B that for strong 
shocks (appropriate to "cold inflow" as in rich clusters 
accreting small clumps and diffuse gas) the approxi- 
mation kT2 ~ — 02/3 + 3A:Ti/2 holds, where (/)2 is the 
gravitational potential energy at r2 — R. 

Inside R, the temperature and density profiles T(r) 
and n{r) are matched to T2 and to n2 by polytropic 
profiles or by their isothermal limit. We numerically 
compute such profiles using the hydrostatic support 
of pressure against gravity (see Appendix B); for defi- 
niteness, we use the Navarro et al. (1996) representa- 
tion for the potential and for the velocity dispersion 
(which varies slowly in the relevant region). 

Let us consider first for reference the simple analyt- 
ical approximation provided by the standard isother- 
mal model 

n{r)/n2 = [p{r)/p2f , (5) 

(Cavaliere & Fusco-Femiano 1976) with the canon- 
ical exponent (3 = fimHC / kT2; here a is the one- 
dimensional velocity dispersion at R, and ^ is the 
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average molecular weight. For the purpose of the an- 
alytical approximation we may use the strong shock 
limit for T2, and find numerical values for (3 ranging 
from about 0.5 for groups to 0.9 for rich clusters (a 
trend consistent with data being collected by M. Gi- 
rardi and coauthors, private communication). This 
implies that Rx /Ris larger in the former than in the 
latter (see CMT97). 

The full numerical computations using the expres- 
sion for T2 given in Appendix B, and cr(r) and (j)(r) 
from Navarro et al. (1996) confirm this trend and 
give results shown in fig. 1. Note that even in the 
isothermal case the; emission-weighted temperature 
(integrated along the line of sight) declines outwards 
but only very slowly. 

The observed stronger decline requires a polytropic 
equilibrium, where the run of T(r) steepens with the 
index 7 increasing from 1. In the Appendix B we 
recall the basic relations, and show that the varia- 
tions induced in the volume-averaged luminosity by 
increasing 7 are small. Thereafter, we adopt the value 
7 = 1.2, which also yields for rich clusters an inte- 
grated baryonic fraction 0.15 out to R. The result 
for the emission-weighted T{r) (see fig. 1) is a mild 
decrease out to r 1 Mpc in agreement with the 
observations (Hughes, Gorenstein & Fabricant 1988; 
Honda et al. 1997; Markevitch et al. 1997), followed 
by a sharper drop as indicated by state-of-the-art sim- 
ulations (e.g., Bryan & Norman 1997). Note that in 
fig. 1 the discontinuity at the shock has been smeared 
out to a smooth drop by the finite resolution, taken 
at 100 kpc for comparison with the simulations and 
with the forthcoming observations. 

2.2. Merging histories and the L-T correla- 
tion 

The luminosity L on v?{r)/n\ p^/^ is sta- 
tistically affected by the merging histories as follows. 
For a cluster or group of a given mass (or temper- 
ature), the effective compression factor squared (5^) 
is obtained upon averaging eq. (4) over the sequence 
of the DM merging events; in such events, T2 is the 
virial temperature of the receiving structure, and T\ 
is the higher between the stellar preheating temper- 
ature and that from "gravitational" heating, i.e., the 
virial value prevailing in the clumps being accreted. 

All that is accounted for in our model using Monte 

Carlo simulations of the hierarchical growth of the 
DM halos; these are based on merging trees corre- 



sponding to the excursion set approach of Bond ct al. 
(1991), consistent with the Press & Schechter (1974) 
statistics (see CMT97). The averaging procedure is 
dominated by the events occurring within the last few 
dynamical times; it results in lowering {g'^) compared 
to g"^, because in many events the accreted gas is at 
a temperature higher than the minimum preheating 
value Tu ~ 0.5 keV. In addition, an intrinsic variance 
is generated, reflecting and amplifying the variance 
intrinsic to the merging histories. 

The net result is shown in fig 2, and commented in 
its caption. In agreement with the observations (Edge 
& Stewart 1991; Mushotzky 1994; Tsuru et al. 1996; 
Ponman et al. 1996) , the shape of the average L—T re- 
lation flattens from L oc at the group scale (where 
the nuclear energy from stellar preheating competes 
with the gravitational energy from infall) to L oc 
at the rich cluster scales. At higher temperatures the 
shape asymptotes to L oc T^, the self-similar scaling 
of pure gravity. Notice the intrinsic scatter due to 
the variance in the dynamical merging histories, but 
amplified by the dependence of L. 

The average normalization formally rises 

like/?i/2(;j)^ 

where p is the effective external mass density which 
increases as (1 -|- z)^ (Cavaliere & Menci 1997) in fil- 
amentary large scale structures hosting most groups 
and clusters (see Ramella, Geller & Huchra 1992). 
This implies a factor 1.3 at 2: = 0.3, consistent with 
the observations by Mushotzky & Scharf (1997). Fur- 
ther weakening of the z-dependence will comes from 
the increasing depth of the central (/)o for distant 
structures of given M, as predicted by Navarro, Prenk 

6 White (1996). 

2.3. The adiabatic models 

At the other extreme, the models by Kaiser (1991), 
and Evrard & Henry (1991) obtain from PE under 
two limits, appropriate only for rapidly expanding 
universes, as we discuss below. The first limit cor- 
respond to no currently active merging, with shocks 
moving outward and vanishing. In such conditions, at 
the boundary T2 « Ti holds; with Ti staying nearly 
constant after the dynamical freeze out, this implies 
g = n2/ni « 1. Correspondingly, the central density 
scales approximately as Uo oc ni (To/Ti)^/''''"^). 

The "adiabatic" models require also a second limit, 
concerning the internal gas distribution. The value 

7 = 5/3 is taken at the center, but an isothermal (3 
profile is assumed (with a fixed /?), based on a King- 
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like DM distribution. A constant baryonic fraction 
at R is then required, and this forces the core radius 
to scale as oc Mi/3-i/3/3. Thus L oc T^+3{2-i/m/2 
obtains, with the normalization p^^^~^^'^{z). 

Finally, the value of (3 is chosen as an input. The 
choice f3 = 2/3 for both clusters and groups leads 
to the model of Evrard & Henry (1991), in which 
L oc T^'^l'^ obtains, with constant normalization. The 
choice /? = 1 leads to the somewhat different model 
of Kaiser (1991), in which L oc T^-^ obtains, with the 
normalization anti-evolving like (1 + z)~^l'^\ this will 
become even more negative when the evolution of 
is taken into account, so conflicting with the data of 
Mushotzky & Scharf (1997). 

2.4. ICP models and cosmology 

Shock and adiabatic models can be characterized 
in terms of entropy (see Bower 1997). Actually, the 
modes of entropy production and distribution corre- 
late with the global dynamics. 

Collapses, merging and the induced shocks are cur- 
rently ongoing in the critical universe, so that strong 
shocks form close to the virial radius. Entropy is con- 
tinuously generated in the outer regions, so that its 
radial distribution is raised outwards. Then the effec- 
tive 7 will be close to one, leading to a roughly flat 
T(r) inside the shock. The density is determined by 
the boundary conditions after eqs. (4). Shocks are 
weaker in groups, the density profile is shallower, and 
the L — T relation steeper. 

Conversely, in an open universe most dynamical ac- 
tion is moved back to early times: merging and mix- 
ing occurred early on, and then subsided; shocks had 
time to expand beyond R and weaken; correspond- 
ingly, the accretion petered out under nearly adia- 
batic conditions for groups and for clusters as well. 
Then the effective 7 is closer to 5/3, and this may be 
used to roughly scale the central densities with differ- 
ent virial temperatures T to obtain L oc T'^ '"' _Rj^. The 
two adiabatic models adopt additional, and different, 
assumptions concerning Vc or Rx, i.e., rc oc T~°-^^, 
or Rx = const, as discussed in §2.3. 

For open cosmologies with ilo « 0.5, or for flat 
ones with Q\ = 1 — Qq, the present deceleration is 
intermediate between the two above cases, and the 
applicability of the shock or of the adiabatic model is 
not so clearcut; we shall consider both, finding similar 
results as is expected. 



3. Statistics of condensations, and X-ray ob- 

servables 

We use the dark mass m, = M/Mo normalized 
to the characteristic mass Mo = 0.610^^ Hq h^^ Mq 
defined in the hierarchical clustering theory, see Ap- 
pendix A and the analytical details given in Appendix 
C. The X-ray temperature T reads: 

T = Tom^/^{l+z) . (6) 

On the basis of §2, for the shock model the bolometric 
luminosity is given by 

L = Lo {g\T)) n\r)lnl {T/Tof {l + z); (7) 

the ratio rfi{r)/n\ (integrated over the cluster vol- 
ume) and the factor {g^) (averaged over the merging 
histories) have been derived in §2.2. 

The constant kTo takes on the value 4.5^0 keV 

(see Appendix C) . The luminosity Lo is calibrated on 
the height of the observed, local L — T correlation 
(see fig. 2) rather than computed a priori, in view 
of the subtleties discussed in §5. At 4.5 keV we find 
Lo = 1.6 (1 ± 0.25) h-^ erg/s. The height of 
the local luminosity function provides an independent 
value for the normalization, which we find consistent 
with the former to within 15%. 

The statistics at different z of DM halos in the HC 
theory is provided by the standard Press & Schechter 
(1974) formula, which in comoving form [po being the 
local cosmological density) reads: 

^ /2 ^cPo , diner , TO~^ ■ ■ 

(8) 

where 5c is the critical threshold for the collapses of 
the density perturbations (depending weakly on the 
cosmological parameters); on the other hand. D{z) 
is the linear growth factor for the density perturba- 
tions, sensitively depending on the cosmological pa- 
rameters (see Appendix A). The linear, time-evolved 
mass variance a{m) D{z) is usually represented in the 
form a%m~°'D{z) where a = (rig + 3)/6 contains the 
effective slope ne{M) w — 1.3-; — 2 of the power spec- 
trum (l^fep) at scales ~ 10-^ 1 Mpc. 

It is characteristic of eq. (9) to comprise two kinds 
of evolution: the number increase oc D~^{z) at the 
\ow-M end, and the shift toward smaller M of the 
upper exponential cutoff. Such dynamical evolutions, 
modulated by cosmology, must combine with the ICP 
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evolution to yield elosely eonstant luminosity func- 
tions as observed. We conservatively adopt the Press 
& Schechter rendition of the hierarchical clustering, 
keeping in mind its problems and limitations (see 
Bond et al. 1991, Cavaliere et al. 1993) and deferring 
discussions to §4.4 and 5. 

From such statistics of dark halos using eq. (8) we 
compute the luminosity function N{L) = N{M)dM/dL, 
the expected flux counts N{> F) of X-ray clusters, 
and their contribution to the soft XRB. The latter 
two observables read (see Appendix C) 



N{>F)=Rh dm dujdi{z)N{m,z)dm 



(9) 



Rh dl{z)N{m) dm [n(r)/ni 



hv 



Lom4/3(l + z)5/2 
47rr(0.6) A;Tom2/3 VkTom^/^ 



„2/3 



(10) 

Here the lower limit m^ is set by the requirement that 
T{m, z) > 0.5 keV, the effective lower bound for group 
temperatures from preheating; smaller masses would 
correspond to galaxies, where the amount of diffuse 
baryons and the emission (per unit total mass) drop 
sharply (Fabbiano 1996). 

The limit mi in cq. (10) is set by the maximum 
between T(m,z) > 0.5 keV and the limiting flux of 
cluster or group sources, for which we conservatively 
adopt 410-^"* crg/s cm . On the other hand, TOi 
also constitutes the upper limit for the the unresolved 
sources contributing to the XRB. 

Thus a complementarity relationship holds between 
the counts and the contribution to the XRB. If one 
limits the number of resolved sources in the counts by 
assuming, e.g., a stricter siuface brightness selection, 
as discussed later on, then all sources pronounced un- 
resolved will contribute to the XRB. The joint con- 
sideration of those two observables is thus expected 
to give robust constraints. 

4. Results 

In the luminosity functions N{L, z) the luminosi- 
ties are reduced to the ROSAT band for comparison 
with the data. The derivation of the temperature 
function N{T,z) involves only eqs. (9) and (7), cor- 
responding to the mere passive role played here by the 
ICR 



The results depend on cosmology, and are sensitive 
to the values of the normalization Lq and of the am- 
plitude fTg. To a weaker extent, they depend also on 
the full shape of the power spectrum, which is deter- 
mined in turn by fio and (weakly) by fi^j but does 
not depend directly on fl\ (although its normalization 
does). 

We note that Lg and as decrease with decreas- 
ing, and by themselves tend to decrease all numbers; 
however, this is delicately balanced by the increase of 
the distances and by the slower (negative) evolution, 
as discussed in more detail in §4.4. With fi;^ ^ the 
amplitude as is larger, and this also enters the results 
as discussed in §4.3. 

As for the cosmological parameters, we conserva- 
tively adopt combinations of h and Oq yielding for the 
present age of the universe = 13 ± 2 Gyr (see Os- 
triker & Steinhardt 1995). On the other hand, many 
X-ray measurements in clusters give a considerable ra- 
tio « 0.15-^0.20 of the baryons to dark matter; as rich 
clusters are likely to constitute fair samples of the uni- 
verse, an abundance ratio fls/^o = 0.05 ±0.02 
is indicated (White et al. 1993; White & Fabian 1995; 
Markevitch et al. 1996). So for ~ 0-3 a sufficient 
value Os = 0.0125 ± 0.0025 h''^ is predicted by the 
standard cosmological nucleosynthesis with canonical 
abundances of light elements (Walker et al. 1991); 
but in the critical case (with h = 0.5) this must be 
stretched up to = 0.15 ^ 0.20. 

The full a{m) for CDM cosmogonies in different 
cosmologies, normalized to the four-year COBE re- 
sults (Gorsky et al. 1996), are given by Bimn & White 
(1996), and White et al. (1996). We focus on three 
popular CDM cosmologies/cosmogonies, which pro- 
vide acceptable values for ag,: the tilted (rip = 0.8) 
spectrum in a critical universe with high baryonic con- 
tent (TCDM); the scale-invariant spectrum at large 
scales {rip = 1) either in a flat imiverse with ilo = 0.3 
and f^A = 0.7 (ACDM), or in an open universe with 
Qo ~ 0.5 (OCDM). In the last subsection, the full 
set of CDM cosmogonies will be discussed in a more 
synthetic way. 

Our results will be compared with the data from 
two recent surveys with ROSAT: the relatively local 
Brightest Cluster Sample by Ebeling et al. 1997, and 
the higher z sample by Rosati et (1997). 
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4.1. Tilted CDM with high baryon content in 

the critical universe 

For the TCDM, we adopt the tilted primordial 
spectral index rip = 0.8 and the amplitude ag = 
0.66(1 ± 0.08), with a high baryonic fraction Ub = 
0.15, in the critical universe with Hubble constant 
h = 0.5. 

The tilt is chosen following White et al. (1996) so 
as to minimize one of the main problems of the Stan- 
dard CDM, namely, the excess of small-scale power, 
still retaining a value for h rather low but still not 
inconsistent with current observations. In addition, 
such a cosmogony includes the high baryonic fraction 
referred to above. We recall that such parameter set 
is hard pressed in terms of the low value of h ~ 0.5, 
and of the primordial abundance implied for light el- 
ements, with that of He exceeding many recent mea- 
surements; the debate is still hot on the related issue 
of the D/H ratio, with recent signs of convergence, 
see Tytler, Fan & Buries (1996); Songaila, Wampler 
& Cowie (1997). 

After §2.4 we consider here only the PE model with 
the full range of shock strengths. In fig. 3a we com- 
pare the local luminosity function with the current 
data; it is seen that an acceptable fit obtains only at 
the lower 2 standard deviations in ag- It is also seen 
that the evolution of N{L, z) is predicted to be virtu- 
ally nil out to z « 0.8, consistent with data by Rosati 
et al. (1997). 

The evolution of the temperature function is char- 
acterized by a fast decrease out to a moderate z in the 
number of high-T (massive) clusters. This is within 
the reach of SAX (Piro et al. 1997) and XMM (Mason 
et al. 1995), and is shown in fig. 3b. 

The counts are shown in fig. 3c. The fit of the pre- 
dicted counts to the bright data reflects the accept- 
ability of the fit to the local luminosity function. Note 
that the slope of predicted counts is sufficiently flat 
to fit both the bright and the faint data by lowering 
as to within the COBE uncertainty. Alternatively, a 
similar result is obtained on using the central value 
of as from COBE, but the higher baryon abundance 
ne/no = 0.20. 

Our computation of the soft XRB (see fig. 3d) 
comprises, as said above, the sources fainter than 
2 20-14 erg/s cm^. Our predictions are compared 
with the residual XRB from ROSAT, once the AGNs 
contribution (rs 70 % at 1 keV, see Hasinger 1996) has 
been subtracted out. The computed curve lies below 



the upper bounds. 

4.2. Open CDM with flo ~ 0.5 

As said, open universes with flo ^ 0.3 apparently 
constitute the simplest way out of the baryonic crisis. 
However, it is seen in fig. 4 that COBE-normalized 
cosmogonies with Slo < 0.4 using the adiabatic model 
for the ICP, suffer from the fatal flaw of severely un- 
derpredicting both the local functions of temperature 
and luminosity. The blame stays mainly with the 
dynamical sectors; spcciflcally, the low spectral am- 
plitude as yields a severe deficit in N{L,0), and cor- 
respondingly a deficit in the bright counts. Similar 
dynamical reasons yield in aimed N-body simulations 
a percentage of complex morphologies lower than ob- 
served in the local clusters (West, Jones & Forman 
1995; Mohr et al. 1995). 

On the other hand, intermediate values of Oq are 
still a possibility. After Liddle et al. (1996), the range 
allowed to the class of open cosmologies by a set of 
observational constraints, including /i > 0.5 and the 
COBE normalization, is narrowed down to vicinity of 
flo ~ 0.5. We focus first on the representative OCDM 
cosmogony, with = 0.5, h = 0.65 with fls = 0.07, 
which yield as = 0.76 (1 ± 0.08). 

Here it is not easy to decide a priori whether the PE 
or the adiabatic model applies better to the ICP, so 
we use both, ending up in similar results as expected. 
These are shown in figs. 5a, 5b, 5c, 5d. It is seen 
that the local distributions are well fitted, but the 
integrated observables show excesses over the data. 

These persist when other values of fig around 0.5 
are used, and when the uncertainties in erg and Lq are 
considered, as discussed in §4.4. 

4.3. CDM with flo = 3 in flat geometry 

"Intermediate" conditions for the cosmic deceler- 
ation also obtain when fi^ ^ is accepted, with a 
flat geometry as in most variants of inflation. The 
values flo = 0.3 and fl\ = 0.7, with fls = 0.05 and 
h = 0.7, match many observational evidences (see 
Ostriker & Steinhardt 1995). Following Klypin, Pri- 
mack and Holtzman (1997), here the normalization is 
as = 1.1(1 ±0.08). 

Here as in the previous case we have to consider 
both the PE and the adiabatic model. We show in 
fig. 6a and 6b the local luminosity and temperature 
functions, while in figs. 6c and 6d we plot the pre- 
dicted counts and the contribution to the soft XRB. 
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4.4. A synthetic presentation 

Here we give results covering the full set of COBE- 
normalized CDM cosmogonies. For a synthetic pre- 
sentation, in figs. 7a and 7b we compare the predic- 
tions for the counts at bright {F = 2.510^^^ crg/s 
cm^) and at faint fluxes {F = 410~^^ erg/s cm^) 
with the data. We also show the effects of varying ag 
within the imcertainty associated with COBE data, 
and Lo within the minimum dispersion intrinsic to the 
L — T relation. Note in the figures that the strips cor- 
responding to the uncertainties narrow down at the 
upper edge, because as (entering inversely the expres- 
sion for AN/N) increases with Clg. Similarly, AN/N 
is larger for the TCDM cosmogony compared with 
Standard CDM, due to the smaller value of erg- 

Our results agree with Liddle et al. (1996) in rul- 
ing out flo < 0.45 and Oq > 0.55. Even in the re- 
maining range our results compare critically with the 
data, because on fine-tuning fio toward « 0.4 the lo- 
cal luminosity function and the bright counts turn out 
to be underestimated; on the other hand, as soon as 
Oo ~ 0.5 is approached, an acceptable fit to the local 
luminosity function is recovered, but an excess in the 
faint counts is generated (see fig. 5), especially with 
the adiabatic model. In summary, agreement with 
both the brigth and the faint data is at best marginal; 
the underlying reason is that in open cosmologies long 
lines of sight and slow dynamical evolution conspire 
to yield a slope of the counts too steep to account for 
both faint and bright counts. 

We note that such a slope would be even increased 
on considering that the formation z is always larger 
than z at observation, which has the affect of steep- 
ening, if only slightly, the luminosity functions (Cava- 
liere, Colafrancesco & Menci 1993; Kitayama & Suto 
1997). Note also that the addition of Qx = 0.7 to 
flo = 0.3 implies a higher value of erg and hence a 
higher levels of faint counts, though lower than in the 
flo « 0.5 case. 

In figs. 7 we also show that the counts in TCDM 
critical and in ACDM cosmogonies/cosmologies can 
be made consistent with the observations on consid- 
ering not only their uncertainties, but also those in the 
present COBE normalization and the intrinsic uncer- 
tainty in the L-T relation, see the discussion following 
eq. (7). 

Overall, a common feature of all the above models 
based on canonical hierarchical clustering, is consti- 
tuted by some excess in the counts; only in the crit- 



ical and in the flat geometry this can be brought to 
consistency with the data. This may indicate some 
non-trivial incompleteness in the canonical hierarchi- 
cal clustering, worth keeping under scrutiny. 

We also show in fig. 8 the results for the contribu- 
tion to the soft XRB, e.g., at S « 1 keV. Once again, 
CDM with 0.55 < flo < 1 is ruled out, while flo ~ 0.5 
is marginal also in this respect. 

Could excess faint counts be reduced by consider- 
ing a stronger selection due to surface brightness? On 
the contrary, we stress that the complementarity with 
the contribution to XRB makes any such excess even 
more significant. In fact, a solution cannot be sought 
in terms of surface brightness selections without in- 
creasing the excess contribution to the XRB from the 
unresolved groups and clusters. 

5. Conclusions and discussion 

In this paper we have computed the X-ray observ- 
ables for groups and clusters of galaxies. As antici- 
pated in the Introduction, use ratlica' than con- 
tinuous and possibly degenerate parametrizations - 
only discrete combinations of physical models appro- 
priate for the Dark Matter and for the Intra-Cluster 
Plasma. 

We first list our results, and then discuss them in 
detail. 

We have developed the punctuated equilibria (PE) 
model for the ICP state and dynamics. This is com- 
prised of the following two components. 

As for single, clusters, we have used a polytropic /?- 
model which yields temperature profile T(r) (see fig. 
1) in good agreement with the observations. We pre- 
dict the ICP density profile n{r) and the brightness 
profile to be flatter for groups than for clusters, cor- 
responding to a larger extension of the ICP relatively 
to their gravitational radii. 

As for statistics, we convolved the ICP equilibria 
with the histories of DM halos, and predicted the L — 
T correlation to take the form shown in fig. 2, in 
agreement with the data. In addition, we predicted 
an intrinsic variance with the minimum value also 
represented in fig. 2. 

Based on our PE model we then proceeded to com- 
pute for various standard cosmological frameworks 
the local and the evolved luminosity functions of 
galaxy clusters, that we compared with the data (fig. 
3a, 4, 5a, 6a). We derived also the number counts 
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(fig. 3c. -50, 6c). the z-distribiitions (fig. 9) and 
the contribution to the soft X-ray background (fig. 
3d, 5d, 6d). Our results are summarized in figs. 7a 
and 7b; these show that the set of acceptable cos- 
mogonies/cosmologies is restricted to three disjoint 
domains: 0.4 < flo < 0.5 for the standard CDM; 
n = 1 for the Tilted CDM: fi,, » 0.3 for CDM in 
flat geometry. In fig. 9 we summarize the confidence 
levels at which the data are matched. 

We next proceed to discuss in detail the results 
listed above. 

5.1. ICP state in evolving DM halos 

The ICP state in the hierarchically evolving grav- 
itational wells constitutes the focus of our new ap- 
proach. We propose that such state follows suit, 
passing through a sequence of punctuated equilibria 
(PE) that we compute semi-analytically. These com- 
putations comprise: the merging histories of the DM 
potential wells, obtained with a large statistics from 
Monte Carlo simulations of the hierarchical cluster- 
ing; the inner hydrostatic equilibrium disposition, up- 
dated after each merging episode; and the boundary 
conditions provided by strong and weak shocks, or 
even by a closely adiabatic compression, depending 
on the ratio of the infall to the thermal energy in the 
preheated external medium. 

The results of our model depend on two param- 
eters, the external temperature Ti and density ni, 
which are not free. Specifically, we use for Ti the 
lower bound Tu — 0.5 keV provided by the litera- 
ture on stellar preheating; in the merging events the 
effective Ti is the virial temperature of the incoming 
clumps, when this is larger than 0.5 keV. The value 
of ni for rich clusters is related to the DM density by 
the universal baryonic fraction fls ~ 0.15. 

Note that our PE model does not require strict 
spherical symmetry, but rather that the residual in- 
ternal velocities be smaller than the inflow velocity. 
So they can include merging episodes ranging from 
nearly isotropic accretion of small clumps and diffuse 
gas, to anisotropic coalescence of comparable clumps 
along fllaments of the large scale structures. 

The expression of the bolometric luminosity is pro- 
portional to = (n2/ni)^, the square of the density 
jump at the bounding shock. The average of such 
factor over the merging histories is what gives to the 
statistical L — T correlation the curved shape shown 
in fig. 2. For rich clusters we obtain L oc T^. This 



flattens to L oc for larger T, corresponding to 
the saturation of the shock compression factor, i.e., 
g{T/Ti) ^ 4 when T > Ti. At the other end, the cor- 
relation steepens to L oc in the group range, where 
T/Ti ~ 1 and the shocks are substantially weakened 
by the preheating temperature in the infalling clumps. 
The amplitude of the L — T correlation rises gently 
proportionally to {1 + z) oc pi{zY^'^ where pi is the 
density in the large scale structures hosting clusters 
and groups. 

In addition, our PE approach predicts an intrin- 
sic variance of dynamical origin due to the different 
merging histories, and built in the factor g'^. Such 
variance constitutes a lower bound, in view of addi- 
tional contributions from the variance in the ambient 
density, and from the central luminosity associated 
with cooling flows (Fabian et al. 1994; White, Jones 
& Forman 1997). 

5.2. Contact with hydro dynamical simulations 
and with observations 

The PE model includes, in a simplified semi-analytical 
form, compression and shocks at the boundary with 
the surrounding environment, which is modulated in 
density by the large scale structure and in tempera- 
ture by the stellar preheating. 

Simulations now clearly show the shocks occurring 
also in major merging events (Schindler & Mueller 
1993; Roettiger, Stone & Mushotzky 1997). The in- 
clusion of the Rankine-Hugoniot conditions rises the 
internal temperature at the expenses of the inflow ve- 
locities. Complex features are found like residual ki- 
netic pressures, and unmixed hot spots in the temper- 
ature distribution; but over times of about 2 Gyr the 
residual kinetic pressure over cluster scales reduces to 
less than 20 % of the thermal one. 

On the other hand, in the hierarchical clustering 
such major events are rare; our Monte Carlo simula- 
tions give a probability <J 20% for very asymmetric 
events with mass ratios 2:1 or larger occurring within 
2 Gyr from the cluster observation. In addition, in 
such events the ICP temperature in the infalling sub- 
cluster is comparable to the virial value in the main 
cluster. Such major events with their low frequency 
and large Ti yield a minor contribution to the statis- 
tical {{n2/niY). 

Our semi-analytical model describes only crudely 
these transient if conspicuous features, to focus on the 
lesser and more symmetric events which contribute 
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the most to the L — T relation. At the extreme 
of spherical accretion of loose gas the simulations 
(see Takizawa & Mineshige 1997) show in detail that 
shocks also form and expand slowly, to leave inside 
a declining temperature profile and a steeper density 
profile. 

Our PE model yields temperature profiles decreas- 
ing as shown in fig. 1. These agree with the published 
data (see Hughes et al. 1988; Honda et al. 1997; 
Markevitch et al. 1997); they agree also with the re- 
sults from state-of-the-art simulations (Bryan & Nor- 
man 1997) obtained by running on supercomputers 
advanced 3D Eulerian codes with adapting mesh and 
reliable shock capturing methods. While the high- 
resolution simulations are limited (for now and for 
some time to come) to the condition of no stellar pre- 
heating suitable only for very rich clusters, our model 
includes the effects of stellar preheating of the outer 
gas over the whole range from groups to clusters. 

5.3. Constraining cosmology 

With the ICP state so described, we proceeded to 
constrain the cosmological parameters. After the ob- 
servations by Rosati et al. (1997), the main rule of 
the game turns out to be as follows: the dynamical 
evolution contained in the standard Press & Selieehter 
formula (eq. 9) must combine with the evolution of 
the L — T correlation and with cosmology to yield 
closely non-evolutionary N{L,z). 

We stress that such combinations are severely se- 
lected in our approach. In fact, strong shocks are 
common in the critical cosmology, where accretion 
and merging activity are currently ongoing, and there 
shocks apply in full; these include also weak shocks 
for small groups with virial temperatures below 1 keV, 
and for those rich clusters which merge with compara- 
ble clumps. Closely adiabatic compressions, instead, 
prevail for all structures in very open universes with 
high formation redshifts, and little or no strong shocks 
and mixing at present. This defines the domain of 
applicability of the two adiabatic models by Kaiser 
(1991) and Evrard & Henry (1991). 

In detail, our results are as follows (see fig. 7). 

We confirm that Standard CDM in the critical cos- 
mology is definitely ruled out on account of its over- 
production of local clusters and of their considerable 
positive evolution. 

We also rule out open cosmologies with fio < 0.3, 
the simplest way to enforce little evolution and also 



to provide a solution to the baryonic crisis. In fact, 
these cosmologies when COBE-normalized yield a se- 
vere deficit in the local luminosity function and in the 
source counts. 

Thus, we investigated more elaborate solutions: the 
critical universe with tilted CDM and a high baryon 
content; and two cases of intermediate deceleration, 
comprising CDM with canonical nucleosynthesis ei- 
ther in mildly open universes with Oq « 0.5, or in flat 
geometry with Oq = 0.3 and fix = 0.7. 

With the first solution we obtain (marginally) ac- 
ceptable results (see figs. 3a-3c) for the local lumi- 
nosity function N{L, 0), for the source counts, and for 
the contribution to the soft XRB, within the present 
uncertainties of the data and within the variance in- 
trinsic to the theory. We note that the cosmologi- 
cal/cosmogonical sectors by themselves may be tested 
independently, based on the fast evolution with z (see 
flg. 3b) of the temperature distribution in the criti- 
cal pointed out by many authors (see Oukbir 
& Blanchard 1992; Henry 1997); this constitutes an 
important program for the satellites SAX and XMM 

As for comologies with intermediate deceleration, 
here neither the PE nor the adiabatic models for the 
ICP are cogently indicated; thus we considered both, 
obtaining generally similar results as expected. For 
open cosmologies in particular, the results are incon- 
sistent with the observations of the local luminos- 
ity function and of the counts, except for the range 
flo = 0.45 -i- 0.55; even there the counts are excessive 
at more than the formal 99% confidence level (see fig. 
9), the excess being larger for the adiabatic models. 
The excess is due to built-in reasons, that is, the rela- 
tive large amplitude ag and the relatively steep shape 
of the counts, as spelled out in §4.4. Instead, in the 
flo + = ^ cosmology a manageable count excess is 
obtained. 

We note that the Kaiser's (1991) variant of the 
adiabatic models does not yield such an excess for 
flo ~ 0.5, due to its normalization decreasing at high 
z. However, the local luminosity function computed 
from this model overestimates the number of bright- 
est clusters, due to the strong dependence L oc T^-^. 
Moreover, the normalization decrease at high z is 
hardly consistent with the data by Mushotzky and 
Scharf (1997). Finally, it yields a fast, negative evolu- 
tion of N{L, z) barely consistent with the data in the 
survey by Rosati et al. (1997) (the deficit is truly fatal 
in the critical or in the flat cosmology). As luminosi- 
ties larger than some 3 10^^ erg sec~^ are little repre- 
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scmtcd in that survey, a strong test for such a negative 
evolution concerns any deficit at bright fluxes in the 
redshift distribution from a large-area survey. So we 
show in fig. 10 the redshift distribution computed 
also for this model. 

We have conservatively chosen to focus on a lim- 
iting flux F = A 10^^^ erg sec~^ for which the sky 
coverage is nearly 100% (P. Rosati, private commu- 
nication), and any incompleteness is out of question. 
Incompleteness due to surface brightness may be rel- 
evant at fainter fluxes, depending on the cluster and 
group profiles. We plan to treat such issue elsewhere, 
but here we point out that in our approach the im- 
pact of any such incompleteness is limited by the com- 
plementarity between counts of resolved sources and 
contribution to the XRB from the rest. 

Our summary is that many combinations of stan- 
dard cosmogonies/cosmologies with ICP models are 
ruled out. A relatively small set of disjoint cos- 
mologies/cosmogonies survive, as shown by fig. 7: 
n = 1 with tilted CDM and high baryonic abundance 
combined with the Punctuated Equilibria, which is 
marginally consistent with the data; CDM in open 
cosmology with 0.5 < fio < 0.55, which is barely 
consistent using the PE, and even less so using the 
adiabatic models; CDM with fio = 0.3 and fix = 0.7, 
which is consistent using either the PE and the adia- 
batic models. 

So cosmological parameters can be constrained on 
the basis of X-ray clusters, but only up to a point; 
for example, the residual uncertainty in the density 
parameter is AOq/Oo > 20%. 

5.4. What next 

To what extent enlarging the data base on X-ray 
clusters will help in further constraining cosmology? 
Here we argue that the variance intrinsic to the hierar- 
chical clustering, and amplified by the ICP emissivity, 
sets an effective limitation. In fact. fig. 7 shows that 
the present Poissonian error bars in the observed faint 
counts are already smaller than the (minimum) in- 
trinsic variance in the predicted ones. Decreasing the 
former with richer, faint surveys will hardly provide 
a sharper insight into cosmology unless one reduces 
both the uncertainty concerning erg and the larger one 
concerning Lg-, in fact, the two enter with compara- 
ble weights eq. (10), since it is seen that Aag/as acts 
roughly as (n^ + 3)ALo/QLo ~ 0.2ALo/Lo. 

To what point are these reductions feasible? On 



the theoretic;al side, the minimum ALg/ Lg of dynam- 
ical origin may be sharpened by Monte Carlo simula- 
tions so extensive as to provide the full scatter distri- 
bution. But then one must tackle also the enhanced 
emissivity produced or signaled by cooling fiows, cor- 
related with higher ambient densities; this we shall 
treat elsewhere (Cavaliere Menci & Tozzi, in prepara- 
tion). 

On the observational side, one needs a large statis- 
tics for the distributions of L and T; this will help 
in deriving narrower L — T correlations for subsam- 
ples categorized in term of mass deposition rates from 
cooling flows, see White, Forman & Jones (1997). 
Such aim calls for spectroscopic measurements of T, 
which are obviously harder than the bolometric L, 
and require SAX or even XMM. However, we stress 
that such efforts will flnd soon a more proper aim than 
constraining Qo- 

This is because soon MAP (Bennett et al. 1997), 
and subsequently PLANCK (Bersanelli et al. 1996), 
will accurately measure on very large and still lin- 
ear scales not only the perturbation power spectrum 
(from which erg is derived), but also directly flo to 
better that 10%; this will supersede constraints set at 
cluster scales gone non-linear. 

Once the cosmological framework has been fixed, 
the study of groups and clusters in X-rays will resume 
what we submit to be its proper course; that is, the 
physics of systems of intermediate com,plexity which 
is comprised of the DM and of the ICP component. 

With the latter fully understood and the scatter in 
the L — T relation assessed, cluster X-raying will fi- 
nally expose the underlying process of non-linear con- 
densation of DM on scales 1 — 10 Mpc. Then any mis- 
match concerning the number counts or N(L, z) will 
be telling of failures either in the CDM spectra or in 
the current representation of cosmogony in terms of 
the Press & Schechter formula. 

As a relevant example, we recall from §4.4 that 
even the acceptable models we computed tend to ex- 
ceed the observed faint counts, and can be brought 
to consistency only at the lower end of the current 
uncertainty concerning erg. On the other hand, the 
corrections to the Press & Schechter formula currently 
discussed yield a larger number of clusters. For exam- 
ple, Jain & Bertschinger (1995), and Gardner, Tozzi 
& Governato (1998) find that the threshold 5c must 
be lowered from the canonical value 1.69 to 1.5, at 
least at z > 1 if not already at ^ = 0; a similar trend 
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obtains considering that the formation redshift is al- 
ways greater than the observation's as discussed in 
§4.4. If MAP wiU provide definite values of Up, trg 
and flo such as to enhance the excess in the faint 
counts, then the Press & Schechter rendition of the 
non-linear cosmogony will have to be reconsidered. 

Preliminary computations of the source counts and 
the contribution to the soft X-ray background have 
been performed by F. Lupini in his Thesis. We 
are indebted with P. Rosati for communicating his 
data prior to publication and for many informations. 
We thank S. De Grandi, C. Gheller, M. Girardi, S. 
Molendi, L. Moscardini, O. Pantano and D. Trevese 
for helpful discussions. We also thank our referee 
M. Henriksen for several helpful comments, and for 
stimulating us to better focus our exposition. Partial 
grants are acknowledged from MURST and ASI. 



A. Standard hierarchical clustering 

We recall that the variance of the perturbation field 
at a scale R = 27r/fc (associated to the mass scale 
M — pAuR^ /i in terms of the average density p) is 
defined by 



a2(i?) = -i j d'k\5u\^W{kR) , 



(Al) 



where W{kR) is a top-hat filter (sec Bond et al. 1991; 
Lacey & Cole 1993). This is given in terms of the 
power spectrum \5k\'^ = k"pT'^{k) which depends on 
the physics of the early Universe (see Peebles 1993; 
Lucchin, Matarrese and MoUerach 1992) and by the 
subsequent microphysics. 

Results from the COBE/DMR experiment give 
Up ~ 1 ± 0.2 and provide the normalization on large- 
angle scales (see Bunn & White 1996 and references 
therein). The transfer function T{k) depends on the 
nature of DM which is the main constituent of the 
perturbations; for CDM, and given values of the pa- 
rameters flo, h and fls, standard formulae are given, 
e.g., by Sugiyama (1995); White et al. (1996). These 
can be recast into the form a{m) = a^mr'^, where 
a = (rie -I- 3)/6, with being the effective, scale de- 
pending, index in the full power spectrum ( 1 5fe P ) . The 
values of the amplitude at 8h~^ Mpc are given, e.g., 
by Bennett et al. (1996) for different cosmologies. 

The scale 8h~^ Mpc also defines a characteristic 



mass 



Mo 



—Po (8 h-^Mpcf 



0.6 10^^ /i-i Mq , 
(A2) 

which we shall use as our unit mass. 

At any given mass scale, the time evolution is de- 
rived from the linear growth of the perturbations D(t) 
(see Peebles 1993), which depends both on the den- 
sity parameter fig and on the parameter fl\ — A/3 
associated to the cosmological constant A. An expres- 
sion valid in all cases with il.o + — I (as predicted 
by most inflationary models) is given by Lupini (1996) 
in terms of the epoch t, and writes 

-1 1/3 



m = 



(A3) 



where to is the present epoch, and a 
For n — 1 the one obtains D{z) = (1 + z)^^. 

The mass distribution N{M, z) of the condensa- 
tions given in eq. (9) has been derived by Press 
& Schechter (1974), and is discussed by Bond et al. 
(1991). Here we stress that it contains only the linear 
mass variance a{'in) and the threshold for non- linear 
collapse 5c for which the canonical value 1.69 is taken. 



B. Postshock conditions, and polytropic equi- 
librium 

The jump conditions at the shocks are based on 
the Rankine-Hugoniot conservations, and may be ob- 
tained from the implicit expressions given, e.g., by 
Landau & Lifshitz (1959). We work out the explicit 
expression of the post-shock temperature T2 for three 
degrees of freedom and for a nearly static post-shock 
condition with V2 « vi, in the form: 



kTo = 



(l + ^/T+7)2 



7 



lo"^ 20(1 



(Bl) 

where e = \bkTi/ Apmnvi. In a "cold inflow" with 
e << 1 the shock is strong, and the expression simpli- 
fies to kT2 ~ pruHvl/i -f 3fcri/2. The fiow velocity 
vi is set by the potential drop across the region of 
nearly free fall, to read vi ~ —lA(f)2/mH where 4>2 is 
the potential at r = i? (see CMT97). Since 1.4// « 1, 
the above equation may be effectively approximated 
by kT2 ~ -cj)/i + 3fcri/2[|. Instead, for e > 1 the 



^We correct here a numerical error in eq. (3) of CMT97, where 
7/8 appeared instead of 3/2. 
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shock is weak and T2 — Ti is recovered as expected. 

When the temperature profile is polytropic with 
T(r) (X n{r)'^~^, eq. (2) is modified to 

1/2 

(B2) 



Lex TV'/' 



T 



[n{r)/n2 



2+(7-l)/2 



The ratio n{r)/n2 is obtained starting from the 
hydrostatic equihbrium dP/ndr = —d(j)/dr with the 
polytropic pressure P(r) = kT2n2 [n{r)/n2^'^ / ^i. This 
yields (see Cavaliere & Fusco Femiano 1978; Sarazin 
1988) the profiles 



T{r) rn{r) 



T2 



n2 



7-1 



1 



7 



■/3[02-0M] , (B3) 



where = <p/fj,mH cr^ is the normalized potential; we 
use for 0(r) and (T(r) the forms given by Navarro, 
Frenk & White (1996). 

Eq. (B2) reduces to eq. (2) of the text in the 
isothermal limit when 7 — > 1 and T = T2. The 
volume-averaged factor in eq. (B2) differs from that 
for the isothermal case by less than 20% in the full 
range 1 < 7 < 5/3. 



C. Luminosity functions and their integrals 

The X-ray emission of clusters of galaxies is due 
to optically thin, thermal bremsstrahlung of the hot 
(T ~ 10"^"^^ K) ICP in equihbrium with the cluster 
potential wells (Cavaliere, Gursky & Tucker 1971; see 
Sarazin 1988 for a review). 

The virial theorem provides T cx M/R. The virial 
radius R can be expressed in terms of the cluster 
mass M and the density p to read R cx [M/ p)^/^, 
which yields T oc M^/'^p^/'^. According to the stan- 
dard hierarchical clustering, the DM density inside 
clusters is proportional to the background's, so that 
p « 200p oc (1 -t- zY\ then one obtains 



T cx M2/3 (1 + z) 



(CI) 



corresponding to eq. (7) in the text. The proportion- 
ality factors in eq. (CI) are given, e.g., by Hjorth, 
Oukbir & van Kampen (1997) to imply T = 4.5 keV 
for a cluster with M — Mco- 

The bremsstrahlung spectrum tiv) at the frequency 



V (in Llio frame of the source) is given by 



r -hv/kT 

£{v)<x / d\n^{r)^=^Gf {hv/kT) , (C2) 



where n is the ICP density, V is the volume of the 
emission region and G/ is the Gaunt factor, which 
may be effectively aproximated with the function 
{hiy/kT)~"-^. From eq. (CI) and (C2) one obtains 
in terms of the observed frequency I'o = J^/(l + z) 



L 



47rr(0.6) 



kTcofn^/^ 
(C3) 



Here L is the bolometric luminosity 

L^Lo [n(r) /mf m^/^ (1 + zf'^ , (C4) 



where {n{r)/niY = /g^ d^ r [n{r) / / R^ , and Lo = 
1.610^'*/i~^ erg/s corresponds to 4.5 keV calibrated 
to the local L — T correlation. For comparison with 
ROSAT data, eq. (C3) is to be integrated over 
the local range /^E = 0.5 2 keV. The resuh- 
ing luminosity in the band AE' reads L^^Eim, z) = 
w{lS.E^ m, z) L{m, z), where the correction factor is 

w{AE, m,z)= / dL;e-^/'=^(™'^) Gf{E/kT{m, . 

(C5) 

The integral number counts are given by 



iV(> F) = dzdViz) / dmN{m,z) 



(C6) 



where dV = Rh dcu d£{z)Dl{z)/{l + z^ is the cos- 
mological volume subtended by the solid angle dw, 
Dl{z) is the luminosity distance, and Rh dl{z) is the 
line-of-sight element depending on fio and The 
lower mass mi is that corresponding (after eq. C4 
and C5) to the lowest luminosity L^e detectable, at 
any z, by a survey with the limiting flux F . 

The complementary contribution to the soft XRB 
of the sources with F' < F is given by the expression 



dV{z) 



dm i^vq 



inDliz) 



(C7) 



With the use of equation (C3) this yields the equa- 
tion 11 of the text, where the limits mo and mi are 
discussed. 
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Fig. 1. Temperature (top panel) and density (bot- 
tom panel) profiles for a single cluster of 10^^ Mq at 
= in polytropic equilibrium; 7=1 (solid line), 1.1 
(dotted) and 1.2 (dashed), see eq. (B3); the dark mat- 
ter potential is taken from Navarro, Frenk & White 
(1996), and critical cosmology with TCDM pertur- 
bation spectrum is used. The shock strength corre- 
sponds to T2/T1 = 10. The temperature is emission- 
weighted along the line of sight, and smoothed with 
a filter width of 100 kpc. 

Fig. 2.— For the PE model and the TCDM cos- 
mogony, we show the average L — T relation (in terms 
of the emission temperatures) with its 2 — a dispersion 
(shaded region). The average and the scatter of the 
compression factor are computed by convolving in 
eq. (4) with the merging histories. The steepening 
at low T is due to the preheating temperature; this 
is uniformly distributed in the interval 0.5 ± 0.3 keV. 
Data from Ponman et al. (1996), solid squares; David 
et al. 1993, open squares. 

Fig. 3. — 3a: The average local luminosity function 

(solid line), and that evolved out to z = 0.7 (dashed 
line) are computed for the shock model and TCDM 
cosmogony. The data are from Ebeling et al. 1997. 
3b: Same as fig. 3a for the local temperature function 
(solid line), and that evolved out to z = 0.5 (dashed 
line). The data are from Henry & Arnaud (1991). 
3c: The predicted source counts are compared with 
data by Rosati et al. (1997; sohd squares) and by 
Piccinotti et al. (1982; solid triangle). Solid line: 
(78 = 0.66; dotted line: cts = 0.55. 
3d: Contribution of the sources with F < 410"^'' 
erg/cm^s to the XRB, compared with the observed 
values (open stars, Hasinger et al. 1997) with a 70% 
contribution from the sources resolved by ROSAT 
subtracted out (solid squares). 

Fig. 4. — The local luminosity function in OCDM 

with ilo = 0.3 for the shock model (solid line) and for 
the adiabatic model (dashed line). Data as in fig. 2a. 

Fig. 5. — Predictions for OCDM cosmogony with 

n,, = 0.5. 

5a: luminosity functions for the shock model with 

the appropriate L T relation: local (solid line) and 

evolved {z = 0.7, dashed line) 

5b: Same as 5a, for the adiabatic model. 

5c, 5d: Same as in figs. 3c, 3d. Both the shock model 



(solid line) and the adiabatic model (long dashed line) 
are shown. 

Fig. 6. — Same as in figs. 5, but for ACDM cos- 
mogony. 

Fig. 7. — Fig. 7a: The predicted faint counts at = 
4 10~^* erg/s cm^ are shown on logarithmic scales for 
the FULL range of flo in Standard CDM cosmogony. 
The solid line is computed for the amplitude ag of 
the perturbation spectrum corresponding to central 
values of the COBE normalization; the strips enclosed 
between the dashed and the dotted lines represent the 
uncertainties in ag added to the intrinsic variance in 
Lo- The predictions from TCDM cosmogony (open 
stars) and from ACDM cosmogony (open circle) are 
also shown. The horizontal lines correspond to the 
upper and lower error bars of the data by Rosati et 
al. (1997). 

Fig. 7b: Same of fig. 7a, but for the bright counts 

at F = 2.5 10"^"'^ crg/s cm^: here the horizontal lines 
correspond to the errors estimated by Piccinotti et al. 
(1982). 

Fig. 8. — The predicted contribution to the XRB 
(at « 1 keV) of unresolved clusters and groups is 
shown for the whole range of f2o in Standard CDM 
cosmogony. The results for TCDM and ACDM cos- 
mogonies are also shown with same symbols as in fig. 
7. The horizontal lines correspond to the data by 
Hasinger et al. (1997), with the 70% contribution of 
resolved sources subtracted out. 

Fig. 9. — This shows the 99% confidence con- 
tours for both the computed local luminosity func- 
tion (solid lines) and the computed number counts 
(dotted fines), in the L44 — as plane (L44 = 
erg/s). The boxes indicate ±1 standard deviations in 
CTs (corresponding to the COBE uncertainty) and in 
L44 (from our Monte Carlo L-T relation). Data as in 
figs. 3-6. 

Fig. 9a: Shock model with TCDM. This cosmogony 
is consistent with the counts within 2 standard devi- 
ations below (Tg- 

Fig 9b) Shock model with ACDM. This shows accept- 
able fits. 

Fig. 9c) Adiabatic model with OCDM (fi,, = 0.5). 
The 99 % contours for the comits are outside the un- 
certainty box for (Tg and L44. 

Fig. 9d) Shock model with OCDM {ilo = 0.5). The 
counts are consistent within 2 standard deviations be- 
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low as- 



Fig. 10. — Fig. 10a: Rcdshift distribution per 
stcradian computed in the OCDM for bright fluxes 
F > 10^^'^ erg cm^^ sec^^, (the corresponding lumi- 
nosity functions arc given in figs. 3, 5, 6). Solid line: 
shock model; dashed line: for the adiabatic model of 
Evrard & Henry (1991); dotted line: adiabatic model 
in the Kaiser (1991) version. 

Fig. 10b: Rcdshift distribution per stcradian com- 
puted in ACDM for fluxes F > 410"" erg cm^^ 
sec~^. Lines as in fig. 10a. 
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